Classical statistical distributions can violate Bell-type inequalities 
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Abstract 

We investigate two-particle phase-space distributions in classical mechanics constructed to be 
the analogs of quantum mechanical angular momentum eigenstates. We obtain the phase-space 
averages of specific observables related to the projection of the particles' angular momentum along 
axes with different orientations, and show that the ensuing correlation function violates Bell's 
inequality. The key to the violation resides in choosing observables impeding the realization of 
the joint measurements whose existence is required in the derivation of the inequalities. This 
situation can have statistical (detection related) or dynamical (interaction related) underpinnings, 
but non-locality does not play any role. 
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I. INTRODUCTION 

Bell's theorem was originally introduced to examine quantitatively the consequences 
of postulating hidden variable distributions on the incompleteness of quantum mechanics 
put forward by Einstein, Podolsky and Rosen 2I (EPR). The core of the theorem takes 
the form of inequalities involving average values of two observables each related to one 
of the two particles. Bell showed that these inequalities must be satisfied by any theory 
containing local variables aiming to complete quantum mechanics in the EPR sense. The 
assumptions leading to Bell's theorem imply the existence of a joint probability distribution 
accounting for the simultaneous existence of incompatible quantum observables 3|, Ij]. Local 
models forbidding the existence of these joint distributions are therefore not bound by Bell's 
;heorem. Indeed local models violating Bell- type inequalities have already been proposed 
sl, y, 0], but these models are mathematical and abstract. In this work, we show that the 
familiar statistical distributions of classical mechanics may lead to a violation of the relevant 
Bell-type inequalities. Our main ingredients will consist first in choosing specific classical 
phase-space ensembles for 2 particles (distributions constructed to be the classical analogs 
of the quantum- mechanical angular momenta eigenstates), and then in choosing detectors 
impeding the existence of the joint probability distribution. We will consider two types of 
settings involving angular momentum measurements, each setting being closely related to a 
well-known quantum mechanical context. The first setting will consist in a classical version 

n 

of the detection loophole [8[: the relevant Bell inequalities will be seen to be violated when 
the sampling is done on a sub-ensemble, defined by the type of detected events, leading 
to averages computed on a partial region of phase-space over which the joint probability 
distribution cannot be defined. Hence the violation has a statistical underpinning - we will 
show there is no violation if the averages are taken on the entire phase-space. The second 
setting will reveal a genuine violation of the Bell inequalities due to dynamical reasons: by 
including in the detection process a local probabilistic interaction between the measured 
particle and the detector inducing a random perturbation that blurs the particles' phase- 
space positions, the derivation of Bell's theorem is effectively blocked, as only correlations 
between ensembles corresponding to a fixed setting of the detectors can be made. This 
example can be seen as a classical version of the quantum measurement of non-commuting 
observables. 
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II. CLASSICAL ENSEMBLES 



We first introduce the classical analogs of the quantum mechanical angular-momentum 
eigenstates to be employed below. The classical distributions of particles can be considered 
either in phase-space or in configuration space; equivalently, one can also consider the distri- 
bution of the angular momenta on the angular momentum sphere. Let us first take a single 
classical particle and assume the modulus J of its angular momentum is fixed. The value of J 
then depends on the position of the particle in the phase-space defined hy VL = {6, (p,pe,P(i>}, 
where 6 and (p refer to the polar and azimuthal angles in spherical coordinates and pg and 
p^ are the conjugate canonical momenta. Let Pz{^) be the distribution in phase-space given 
by 

Pz,ie,<p,pg,p^) = NS{JM - JzMJ'm - Jo)- (1) 

defines a distribution in which every particle has an angular momentum with the same 
magnitude, namely Jo, and the same projection on the z axis Jz^. Hence pz^ can be consid- 
ered as a classical analog of the quantum mechanical density matrix |jm) {jm\ since just like 
a quantum measurement of the magnitude j and of the z axis projection m of the angular 
momentum in such a state will invariably yield the eigenvalues of the operators and Jz, 
the classical measurement of these quantities when the phase-space distribution is known to 
be pz will give Jq and J^p (see Appendix A). Eq. ([1]) can be integrated over the conjugate 
momenta to yield the configuration space distribution 

p{e, <P) = N sin(^) ^Ji-Jl/sm\e) (2) 

where we have used the defining relations J^{^) = Pif, and J'^{Q) = pi + p^/ sin^6'. Further 
integrating over 9 and and requiring the phase-space integration of p to be unity allows 
to set the normalization constant = Jq/2t{'^. 

There is of course nothing special about the z axis and we can define a distribution by 
fixing the projection Ja of the angular momentum on an arbitrary axis a to be constant 
(in this paper we will take all the axes to lie in the zy plane). Computing the distribution 
Pao = ^{Ja — Jao)^{J ~ Jq) IS tantamount to rotating the coordinates towards the a axis in 
Eq. ([2]). Fig. 1 shows examples of configuration space particle distributions and gives for one 
plot the corresponding quantum mechanical angular-momentum eigenstate (the similarity 
is not accidental, as Eq. ([2]) is essentially the amplitude of the spherical harmonic in the 




FIG. 1: Normalized angular distribution for a single particle in configuration space, (a) Quantum 
distribution (spherical harmonic \YjM{d,(p)\'^)- (b) Classical distribution p^^^Ojcj)) of Eq. 
(c) Classical distribution corresponding to a fixed value of Ja (here 9a = 7r/4). The angular 
momentum and the projection on the z [(a)-(b)] or a [(c)] axis is the same for the 3 plots {J/rj = 40, 
with r] = h, and M/J = 5/8). 

semiclassical regime, see Appendix A). We can also determine the average projection Ja on 
the a axis for a distribution of the type ([2]) corresponding to a well defined value of J^: 



where 6a is the angle (^TS) and the projection of the component of Ja on the y axis vanishes 
given the axial symmetry of the distribution. ^ 

The original derivation of the inequalities by Bell [1] involved the measurement of the 
angular momentum of 2 spin-1/2 particles along different axes. Here we will consider the 
fragmentation of an initial particle with a total angular momentum Jt = into 2 particles 
carrying angular momenta Ji and J2 (we will assume to be dealing with orbital angular 
momenta). Conservation of the total angular momentum imposes Ji = — J2 and = J2 = 
J. Quantum mechanically, this situation would correspond to the system being in the singlet 
state arising from the composition of the angular momenta [jx = 0, mi = —m2)- Classically 
the system is represented by the 2-particle phase space distribution 



where N is again a normalization constant. On the angular momentum sphere the distri- 
bution (jl]) corresponds to Ji and J2 being uniformly distributed on the sphere but pointing 
in opposite directions. This distribution will now be employed for determining averages of 
observables related to the angular momenta of the two particles. 



p^cos9a6{J^{Q) - Jzo)dQ = Jzo COS 9, 



(3) 




(4) 
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III. STATISTICAL VIOLATION OF THE BELL INEQUALITIES 



Let us assume two types of detectors yielding outcomes related to the angular momenta 
of the particles. The first type gives a 'sharp' (S) measurement of Jia only if Jia is an 
integer multiple of some elementary gauge 77, and gives elsewhere. This detection can be 
represented by the phase-space quantity 



where flik are the parts of phase space where Jia = kr] compatible with a detection (see Fig. 
2(a)). The second detector gives a 'direct' (D) measurement of J2b (the projection of J2 on 
an axis b). The corresponding phase-space function is 



In classical mechanics there is no natural unit for quantities having the dimension of an 
action, so J and 77 can be expressed in terms of arbitrary units, and any physical result will 
depend only on the ratio J/rj. We will assume for definiteness that 1] is chosen so that the 
extremal values ±J can be reached. J/t] must hence be either an integer or a half-integer, 
the extremal values in dimensionless units being given by ±L = iJ/r^. For example if 
rj = 2J, the measurement can only yield the extremal values L = ±1/2 {t] = J allows to 
measure ±L = ±1 and 0, r/ = 2J/3 allows ±L = ±3/2 and ±{L - 1) = ±1/2 etc.). Note 
that the particle label 1 or 2 can be attached to the detectors: indeed, we will call '1' the 
particle detected by S and '2' the particle detected by D. 

The classical average E{a, h) = {SaDb) for joint measurements over the statistical distri- 
bution p can be computed from 



with Eqs. (jl]), (jSD and Given the characteristics (jSD of the S detection, Eq. ([7]) is 
actually a discrete sum over the parts of phase-space flik leading to the detection of kr]; this 
can be written by including a delta function under the integral. Eq. (jlj) imposes 6*2 = vr — 6'i 
and 02 = TT ± 01, and Eq. ([7]) becomes 



Sa{fli) = Jiai^i) if € ^ik, Sai^i) = elsewhere. 



(5) 



^^(^2) = J2 • b. 



(6) 






k=-L 
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where we have chosen the z axis to coincide with a to take advantage of the axial symmetry 
imposed by Sa (here the hmiting procedure in the delta function is understated). The \ 
prefactor is the only nontrivial normalisation factor (coming from the integration over Q\). 
We obtain the average as 

E{a, b) = -liL + 1)(2L + 1) cosiOh - 9a), (9) 

which as expected depends solely on the ratio J/rj = L. 



The correlation 



unction employed in Bell's inequality can be obtained in the standard 



or CHSH) form We choose 4 axes a, 6, a', h' (we can assume an S detector is placed 

along a and a', and a D detector along h and h') and determine the average values for each 
of the 4 possible combinations involving an S and a D detector. The correlation function C 
relating the average values obtained for different orientation of the detectors' axes is 

C{a,h,a',h') = {\E{a,b) - E{a,b')\ + \E{a' ,b) + E{a' ,b')\) (L)-^ (10) 

where we have divided by to obtain the CHSH function in the standard form characterized 
by values bounded by ±1. Here the detected values obey the conditions \S/L\ < 1 and 
\D/L\ < 1, so that the usual derivation of the Bell inequalities would lead to 

C{a,b,a',b') <2. (11) 

By replacing Eq. ([9]) in Eq. ( ITOi) . it can be seen that for L = |, 1 and |, there are several 
choices of the axes that lead to C{a,b,a' ,b') > 2. The maximal value of the correlation 
function corresponds to C(0, j, |, ^) = A\/2 and 2\/2 for L = ^ and 1 respectively |l6| . 

The violation of the Bell inequality is due to the fact that we are only including in the 
statistics the measurements for which both the S and the D detectors click. But when an S- 
measurement is made along the two different orientations a and a' that enter the correlation 
function, different and mutually exclusive parts of phase-space are selected, so that the 
different events 

{Sia, D2b}, {Sia', D2b}, {Sia, Afc'}, and {Sia', D2b'} (12) 

are not supported by a common phase-space distribution. As a consequence the quantity 

J Sain^)Dbin2)Sa'in^)Dbin2)pin^n2)dn^dn2 (13) 
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FIG. 2: Setups for the first (a) and second (b) examples investigated in this work. In (a) an S 
detector is placed along the a axis and a D detector along b. The angular momenta, originally 
distributed on the sphere, are constrained to move on the rings (red dotted) corresponding to fixed 
values of the projection on a. (b) shows the (z7y) plane of the angular momentum sphere for 
the L = 1 case (hence J = 3/2); the 3 zones correspond to the projections of the spherical zones 



P'c 



-1 „o 



p[J and If Ja is measured, the presence of J in one of these zones yields the outcomes 



Ja = —1, 0, 1 with unit probability. If J5 is measured and J € any of the outcomes Jb = —1, 0, 1 
can be obtained with probabilities depending on the distribution p^. 



describing the average of simultaneous measurements along the 4 axes becomes undefined. 
However, as we mentioned above, the existence of the joint probability distribution in the 
integrand of Eq. (1131) . or equivalently Uj, of a common distribution for the events fll2p is a 
necessary ingredient in the derivation of Bell's theorem, thereby explaining the violation of 
the inequalities. It is noteworthy that if one includes the entire phase-space in the average 
([7]) instead of the parts of phase-space corresponding to the double-click events, then Eq. 
f lT5]) becomes well-defined. It can then be shown that E{a,b) and C{a,b,a',b') should be 
multiplied by the fraction of phase-space yielding the double click measurements IT]] : as a 
result Bell's inequality would not be violated. From the standpoint of classical mechanics, the 
objection regarding the necessity of including the entire phase-space makes sense, since one 
can envisage in principle a particle analyzer able to detect the particles that have not been 
included in the double-click statistics. The quantum analog of this problem is the well-known 

nn 

detection loophole, pending on the experimental tests of Bell's inequalities [81. 113I|. 
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IV. DYNAMICALLY INDUCED VIOLATION OF THE BELL INEQUALITIES 



Our second setting goes further into the violation of Bell's inequalities by postulating 
a model involving a local probabilistic interaction during the measurement between the 
detector and the particle being measured: we then obtain a violation of the inequality for 
the entire ensemble of particles. Let us take two identical detectors Ti and T2 that give as 
only output the integer or half-integer values k = L, L — 1, ... — L of the projection Jia and 
J26 of the angular momenta of the particles. We choose here L = J/rj — 1/2, from which 
it follows that the maximal readout L is smaller than J; for notational simplicity we put 
7] = 1 (so J, rather than J/t] takes integer or half integer values). We further assume that 
there is an interaction between Ti and particle 1 (and between T2 and particle 2) affecting 
the angular momentum of the particle so that the transition Jia —>■ k is a physical process 
due to the measurement. 

We impose the following constraints on this process (which only involves a single particle 
and its measuring apparatus, hence we drop the indices labeling the particles). 

1. There are distributions such that if J G 

P^ri3ep':) = 6,k'. (14) 

This means that if is measured and we obtain k then we know that previous to the 
measurement J & Pa with unit probability. 

2. Let {Jb)pk be the phase-space average of Jb over the distribution p^, where the direc- 
tions b and a are assumed to be different. If is measured and J G p^, any outcome 
k' can be obtained with a non-vanishing probability P^f'{J G p^). Our main assump- 
tion is that averaging over Tb gives the phase-space average of Jb, i.e. the interaction 
vanishes on average. This constraint takes the form 

L 

(^'')pS= E ^'^'^J e P^) = ( J.),. . (15) 

k'=-L 

Eq. (ITSll also holds if 6 = a but then only Pj" = 1 is non-vanishing hence 

(Ta),. = {Ja)p. = k. (16) 
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We will not be interested here in putting forward specific models of the interaction yielding 
such probabilities; it will suffice for our purpose that a set of numbers Pk verifying Eq. f|T5|) 
and obeying Ylik Pk = ^ can be obtained. We need to specify however the distributions 
obeying Eq. f[T^ . It is convenient to specify in terms of the distribution of J on the 
angular momentum sphere: it can then easily be seen that Eq. (1161) is realized if is taken 
to be the ring centered on the a axis and bounded by k — 1/2 < J a < k + 1/2 (see Fig. 
2(b)). Then a measurement of will yield the outcome k with unit probability: 

Ta = kiik-l/2<Ja<k + 1/2. (17) 

One can of course envisage a distribution p obtained by combining the elementary ensembles 
p\. In particular the uniform distribution on the sphere is the sum of the 2L + 1 spherical 
rings p\, 

and therefore if Ta is measured the probability of finding a given value /c is P =1/ (2L + 1). 
Inversely the obtention of the given outcome k is correlated with J G previous to the 
measurement. With p\ defined in this way [Eq. f|T7j) ]. {Jh) pk is computed straightforwardly 
and Eq. (fT5|) becomes 

(Tfc)^, = fccos(^^,-e„); (19) 

we see again that for correlations involving averages, the knowledge of the individual proba- 
bilities Pj** is not necessary. Note however that for the particular case J = 1 (i.e., L = 1/2) 
the constraints (fT5|) - (fT7|) as well as the normalization of the probabilities impose the values 
of the P^** irrespective of any precise physical process: indeed k can only take the values 
± — 1 /2 from which it follows that 

Pi" = ^ ± ( J,) . = l±kcos {9, - 9,) . (20) 

Let us now go back to the 2-particle problem, assuming the initial phase-space density p 
given by Eq. (jl]). The expectation value E{a,b) = {TiaT2b) is computed from the general 
formula 

L 

E{a,b)= kk'P{T2b = k'nT^, = k) (21) 
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where k and k' run on the possible outcomes. The probabihties of obtaining Tia = k and 
T2b = k' are obtained in the following way. Using 

P{T2b = k'n Tia = k) = P{Tia = k)P{T2b = k'\Tia = k) (22) 

we first determine P{Tia = k) by remarking that the initial distribution p corresponds to 
Ji being uniformly distributed on the sphere. According to the results of the preceding 
paragraph, with the sphere being cut into 2L + 1 equiprobable zones [see Eq. (ITSl)]. we 
have P{Tia = k) = 1/(2L + 1). We also know that an outcome Tia = k corresponds to 
Ji G Pa [Eq. (fl^ ]. From the conservation of the total angular momentum, we infer that 
particle 2 must lie in the zone p~'' defined by /c — 1/2 < —J2a < k + 1/2 [Eq. f|T7|) ]: indeed 
if were to be measured we would be assured of finding T2a = —Tia = —k. Hence the 
conditional probability appearing in Eq. fl2^ is given by 

P{T2b = k'\Tia = k) = P^nJ2 e p~') (23) 

where Pj*" was defined in Eq. ( ITSi) . The sum over k' in Eq. ( !2T1) thus verifies Eq. ( fT5l) and 
having in mind Eq. (JT9l) . the expectation value becomes 

E{a, b)=J2 ^lTi " = ~^^3^^^ ^o^^'' - ^'^)- (24) 

k=-L 

The correlation function is again given by Eq. ( fTOll . since the maximum value detected by a 
T measurement is L, not J. The result given by Eq. (1241) is familiar from quantum mechanics 
- it violates Bell's inequality for L = 1/2 with a maximal violation for C(0, ^, |, ^) = 2^/2. 
As noted for the single particle case, the derivation of E{a, b) does not depend in any way on 
the individual values of the probabilities P^,' but only on the condition (]I5|) regarding the 
particle-measurement interaction. Note that by Bayes' theorem, it is of course equivalent to 
compute P(T2b = k' n Tia = k) from P{T2b = k')P{Tia = k\T2b = k'), ie by assuming that 
T2b = k' is known first. 

The violation of the Bell inequalities is due to the conjunction of two ingredients. The first, 
represented by the constraints (fT4l) - (fT6ll . is relative to a single particle and its interaction 
with the measurement apparatus. The second is the conservation of the angular momentum 
on average. Interestingly the first ingredient is the one that contradicts the assumptions 
made in the derivation of Bell's theorem. The reason is that Eqs. ( |T^ -( fT6l) are incompatible 
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with the introduction of elementary probabihty functions p^'' (Q) such that 

plnj^p'a) = jpi^mpimd^: (25) 

indeed, such probabihty functions would need to depend on the ensemble, giving rise to 
functions of the type pl''{^', Pa)- This is shown for the case L = 1/2 in Appendix B. With 
this point in mind, one can expand Eq. (I2T1) (with Eqs. (1221) . ( ITSl) and (1231) ) as 

E{a,b) = J J2kpl''{^i;pi^)pd^i)dni J B{n2,k)p-\n2)dn2 (26) 

k 

with 

B{Q„k)^J2^'P^'^^^-^Pa')- (27) 

k' 

The dependence of i? on A; is the crucial property allowing to violate Bell's inequality 
(whereas the dependence of p(^^2) on k in Eq. (126!) by itself can be absorbed in the ini- 
tial correlation 5(Ji + J2) provided k = k{Qi)). The dependence of 5 on has nothing to 
do with non-locality or action at a distance. It is a simple consequence of the logical in- 
ference characterizing the conditional probability ( l22l) given the characteristics of the single 
particle interaction with the measuring apparatus, namely the fact that the model allows 
only specific types of correlations: in the single particle problem one can only correlate a 
given outcome with a specific distribution - this happens when the distribution is symmet- 
ric relative to the detector's axis [Eq. (fT4l) ]: in the two particle problem the single particle 
property just mentioned makes only possible the correlation of J2 as a function of Ji in 
terms of the ensembles to which they belong, not in terms of their individual positions. This 
is consistent with the fact that the knowledge of the individual position of J is meaningless 
to compute the observed probabilities, as even the elementary probabilities must depend on 



the ensemble to which the angular momentum belongs [18|. 

Note finally that would B not depend on k (and the elementary probabilities on the 
ensembles), Eq. (!26ll would turn into 

E^^{a,b) = j A{ni)B{n2)p{^i,^2)dnidn2 (28) 

where 
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thereby yielding the famihar form taken by the expectation value in the derivation of Bell's 



theorem. In the deterministic case considered by Bell [10| the functions p^" and p^? are 
either or 1 depending on the individual position of Ji (resp. J2). This implies that 
k = k{Qi), ie a given outcome depends on the position of Ji on the angular momentum 
sphere, and p^^{Q2) = p(^2|Ji) does not depend on A; or a but on Ji = — J2 (hence 
the inclusion of the term 6{3i + J2) in the definition of p{Qi,Q2))- Conversely one may 
assume p{il2\k) = p{il2\^i) in Eq- fl26l) with Pa(^i) Pb i^'^) being probability functions 
different from or 1; then A and B defined in Eq. ( |29|) are not the observed outcomes 
but their averages, and E^^{a,b) is the expectation corresponding to the stochastic case 
considered by Bell. Bell's stochastic case correlates the individual positions of Ji and J2 to 
possible outcomes with definite probabilities. In the present model the random interaction 
forbids to make the correspondence between a given position of the angular momenta and 
a definite outcome; instead the correspondence is between a definite outcome and a given 
ensemble describing the positions of the angular momenta compatible with the outcome 
(of course if the former correspondence is satisfied, so is the latter, but the converse is not 
true). In the latter case, the structure of the expectation value ( l26l) does not allow to 
define a term of the type given by Eq. f|T3l) whereby a single distribution can account for 
several simultaneous joint measurements. It appears indeed that the ensemble dependency 
exhibited by the present model is a necessary feature in order to produce non-commuting 
measurements 1^. In this sense the present model can be seen as a classical analogue of 
the quantum measurement of two non-commuting observables (such as Jia and Jia') applied 
to correlations between two particles as originally considered by EPR ^]. 

V. CONCLUSION 

The present results show that averages obtained with 2-particle classical distributions 
constructed to be the analogs of quantum mechanical eigenstates can violate Bell's inequal- 
ities. The violation does not involve nonlocality but statistical or dynamical processes that 
impede the existence of joint probability distributions or the correlation between individual 
values of the variables as required by Bell's theorem. Possible implications on the role of the 



Be 



1-CHSH argument as a marker of quantum nonlocality, which has recently been criticized 



14l |. will be examined elsewhere 12 1. 
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APPENDIX A 



The scheme we are employing to contruct the classical distributions rests on the well- 
known analogy between the classical Poisson brackets and the quantum commutation rela- 
tions in the density matrix formalism. Let G be an operator and \ipg) an eigenstate with 
eigenvalue g. Then the pure-state density matrix pg = \ipg) {ipg\ verifies [pg,G] = and 
Gpg = gpg. In classical mechanics the Poisson bracket of two phase space quantities u{q,p) 
and v{q,p) is a canonical invariant defined by [jj 



dudv dudv , , 

dq dp dp dq ' 

Let p{q,p) be the phase-space distribution and G{q,p) be a function such that {p, G} = 0. 
This means that p is invariant relative to the canonical tranformation generated by G, ie 

{p,G}6Qg = Sp = 0, (A2) 

where Qg is canonically conjugate to G, which is a constant of the motion. Then every 
point of the distribution p will be characterized by the constant value taken by G, denoted 
g. If this is the only constraint imposed on the distribution, p{q,p) will take the form (up 
to a normalization constant) 

p{q,p)=5iGiq,p)-g). (A3) 

In configuration space, the distribution p(g) is obtained by integrating over the values of the 
momentum compatible with a given q, 

p{q) = f p{q,P)dp = I ^^^~^'^ dp (A4) 



dG 




dp 


Pi 



where Pi is the root (assumed to be unique, else a sum is in order) of the argument of the 
delta function. Integrating yields 



dp 

P^'^ = dG 



(A5) 



dqdG ^ 

where S{q, G) is the classical action. The configuration space density is therefore the am- 
plitude of the quantum density matrix element {q\ Pg \q) in the semiclassical approximation. 
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APPENDIX B 



We show that the detection model for a single particle given in Sec. 4 is inconsistent with 
probability functions defined by Eq. fl2S]) in the L = 1/2 case (the one violating the Bell 
inequalities). Take Eq. (12^ with a = b and k,k' = +1/2, 

pj^(j ^pt) = l pi'mptmdn = 1. (Bl) 

Particularizing the general formula (ITTl) to the case L = 1/2, p'^ is the positive hemisphere 
of the unit sphere (since J = 1) centered on the b axis. The result on the right handside 
follows from Eq. (1201) . Eq. (IBip implies that p^{^) = 1 ioi 3 E p^ and consequently 
p^i^l) = 0. Conversely since P+''{J G p^) = 0, we must have p+ (fi) = and p^lQ) = 1 
when J G p^. Now assume that the distribution is instead p+ with a different from the b 
axis. Then according to our model [Eq. (!20|) ] we should have 

p?(j G p+) = 1 p^H^])p+(^])rf^] = cos^ (B2) 

Noting that p^, the positive hemisphere centered on a, is actually composed of two parts. 
Pa ^ Pb ^^'^ Pa ^ Pb write 

p?(J ept)= [ pi{n)pU^)dn + [ pl^iQ)pt{n)dn. (B3) 

But we have seen that = 1 for J G p^ and p+ (f^) = for J G p^, hence 

p?{3 ^pt)= I Ptmn = 1 - (B4) 

/n+nn+ 

which contradicts Eq. flB2p . Hence probability functions obeying Eq. fl25l) do not exist, and 
Eq. (125!) should be replaced by 

Pl^J ^Pa)=l plK^;Pa)plmd^ (B5) 

where the notation p^?(f2; p^) denotes the dependence of the elementary probabilities on the 
distribution. Note also that Eq. ( |25l) does hold if one drops the requirement that p^?(r2) 
should represent an elementary probability: for example the functions p^{^) = Jf, + 1/2 or 
p^{^l) = 2JbH{Jb) fulfill Eq. ( ]B2I) without depending on the distribution, though none of 
these functions is contained in the interval [0, 1] and are thus not probability functions. We 
stress that these features, which put strong constraints on the type of admissible physical 
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models that one could envisage, are relevant to a single particle and its interaction with the 
measurement apparatus. 
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The reader familiar with the Bell inequalities for the quantum measurement of Jia and J26 will 
recognize the similarity of Eq. ([9]) with the quantum expectation value; the only difference 
is that the quantum expectation value is normalized respective to the number of possible 
outcomes (2L + 1) whereas here the normalization is relative to classical phase-space (namely 
the length 2L of the measurement axis) . 
[17] Here this part of phase-space is infinitesimal, since for the sake of mathematical simplicity 
we have modeled the S-detection by a delta function. If we replace the delta functions on 
the angular momentum sphere by narrow rings and spherical caps having a finite surface, the 
fraction of phase-space leading to double-click events becomes finite, and the reasoning as well 
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as the conclusions reached with the delta function modeling hold (although the computations 
need to be made numerically) [l^ . 
[18] It would be of course extremely valuable to understand what kind of physical processes are 
compatible with this type of behaviour (for example the value of the angular momentum in 
this case could represent some time average of an underlying stochastic process, or a space 
average of a field-like quantity distributed all over the ensemble). 
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